Dynamical magnetoelectric effects induced by the Dzyaloshinskii- 
Moriya interaction in multiferroics 
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Abstract. - We study the dynamical interplay between ferroelectricity and magnetism in a mul- 
tiferroic with a helical magnetic order. We show that the dynamical exchange-striction induces a 
biquadratic interaction between the spins and transverse phonons resulting in quantum fluctua- 
tions of the spontaneous ferroelectric polarization P in the ferroelectric phase. The hybridization 
between the spin wave and the fluctuation of the electric polarization leads to low-lying transverse 
phonon modes. Those are perpendicular to P and to the helical spins at small wave vector but 
then turn parallel to P at a wave vector close to the magnetic modulation vector. For helical mag- 
netic structure, the spin chirality which determines the direction of P, also possesses a long-range 
order. Due to the dynamical Dzyaloshiskii-Moriya interaction, the spin-chirality is strongly cou- 
pled to the spin fluctuation which implies an on-site inversion of the spin-chirality in the ordered 
spin- 1/2 system and results in a finite scattering intensity of polarized neutrons from a cycloidal 
helimagnet. 



Introduction. — Multiferroic compounds in which 
the electric and the magnetic order coexist are in the fo- 
cus of current research. Of a particular interest is the 
possibility of controlling the direction of the spontaneous 
electric polarization by a magnetic field [1]. Beside this 
, technological relevance of such a strong interplay between 
■ the magnetic and electric order parameters it is also fun- 
, damentally interesting to understand how such a coupling 
' comes about and what is the microscopic mechanism be- 
' hind the magnetoelectric (ME) coupling in multiferroics. 

Among the family of multiferroics, two different type of 
manganites, RMnOg (R= Gd, Tb, Dy) [2] and RMnaOs 
{R= rare earth, Tb,Y, Bi) [3] play a special role as they 
exhibit different microscopic ME coupling mechanisms: 
In the perovskite multiferroic RMnOa, it's shown experi- 
mentally that the onset of helical magnetic order induces 
spontaneous ferroelectric (FE) polarization, which can be 
well described by the so-called spin-current model [4]. In 
these compounds, the spin-orbit coupling within the d{p)- 
orbitals of magnetic (oxygen) ions produces an electric po- 
larization of the form [5] , P x Sj . Non-collinearity 
of the spins Sj (at sites j) is strictly required in the spin- 



current model. In the main FE phase of RMn2 05 the 
electric dipole moments are directed along the 6-axis, the 
spins however are almost collinear in the ab plane which 
indicates that a spin-orbit-driven mechanism can not be 
the primary source for FE in RMn2 05. As an alter- 
native explanation, the (super)exchange-striction [6, 7] is 
believed to be the origin of ferroelectricity in RMn205, 
P ~ Si • Sj. On the other hand, inspecting carefully the 
dynamical properties of the multiferroics, we find both, 
the antisymmetric Dzyaloshiskii-Moriya (DM) interaction 
and the symmetric magnetostriction play an essential role 
and need to be taken into account. 

Based on the spin-current model, the dynamical prop- 
erties of DM interaction were studied in Ref. [8-10]. A 
collective ME excited mode, so-called electromagnon, was 
theoretically predicted. Moreover, it was found [8] that 
this new low-lying mode is perpendicular to both the 
spontaneous polarization and the helical wave vector. It 
corresponds to a rotation of the spin plane with respect 
to the axis of the helical wave vector, the rotation fre- 
quency is V SJD, where S is the spin value, J is the 
exchange coupling, and D is the magnetic anisotropy. 
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Electromagnons have been detected in RMnOa [11] and 
Euo.75Yo.25MnO,j [12]. seemingly consistent with the theo- 
retical analysis. However, a detailed study of the terahertz 
spectrum of Eui-^^Y^^MnOs [13] revealed that infrared- 
absorption along the spontaneous polarization direction is 
also observable, which is not explained by theory. 

In this paper we show that the dynamical exchange 
striction intrinsically generates bi-quadratic coupling be- 
tween the spin and the transverse acoustic(TA) phonon, 
~ (u^ — Uj-)'^{Si ■ Sj) where is a transverse displace- 
ment at site j . This dynamical coupling does not con- 
tribute any additional static electric polarization but in- 
duces the fluctuation of the electric dipole moment due to 
the low frequency excitation modes of TA phonon. One 
thus has a mode mixing behavior and the polarization 
correlation function follows the soft magnetic behavior of 
the system parallel to the uniform electric polarization P. 
Moreover, in 5 = 1/2 multiferroics, the spin- fluctuation is 
accompanied with an inversion of the local spin and con- 
sequently with an inversion of the on-site electric dipole 
moment according to the spin-current model. The hy- 
bridization between phonons and spins results also in a 
large quantum fluctuation of the spin-chirality which al- 
lows for a finite differential scattering intensity of polarized 
neutrons from a cycloidal magnet LiCu202 [14]. 

Dynamical exchange-striction. — We consider a 

one-dimensional spin chain along the z-axis with a frus- 
trated spin interaction. When the temperature is lowered 
a spiral magnetic structure is realized [14,15]. For such 
a helically ordered magnetic phase the spin-current model 
predicts a uniform electric polarization perpendicular to 
the spin chain. The macroscopic electrical polarization 
P is induced by the condensation of the transverse op- 
tical(TO) phonons, P = — euo- An effective model can 
be introduced to describe the spin-phonon coupling [8] as 
follows 

H = Hs + HjjM + Hp (1) 
H,= Ji{ri-rj)Si-Sj 

+ ^ J^in - rm)Sl ■ Sm 
{lm)„n„ 

Hdm = -A^Ui • [e^ X {Si X Si+i)] 

i 

i i 

where the notation {ij)nn indicates nearest-neighboring 
(nn) i and j, and {lm)nnn corresponds the next-nearest- 
neighboring (nnn) I and m. The competition between the 
nn ferromagnetic interaction (Ji < 0) and the nnn anti- 
ferromagnetic interaction ( J2 > 0) leads to magnetic frus- 
tration and to a spiral spin ordering with the wave vector 
cosQ = — J1/4J2 [15]. Hp is an optical phonon model. 



The spin-phonon interaction Hum originates from a spin- 
orbital (DM) coupling and breaks the inversion symme- 
try along the chain. Minimizing the energy yields the 
condition of the atomic displacement and the local spin- 
configuration 

Ui = ^e^ X (Si X Sj+i) (2) 

Particularly, if the zx helical spins along the chain, i.e. 
Si = ^(siniQ, 0, cosiQ), Eq.(2) leads to a macroscopic 
uniform lattice displacements along the x direction ug = 

sin Qix- 

In the helical spin-ordering phase can't be softened 

through the hybridization between the TO phonons and 
the magnons because k/M ^ JS. The spontaneous FE 
polarization Pj, is frozen at — eug- The fluctuation SP^ 
can therefore be neglected [8]. However, accounting for 
the superexchange striction, TA phonon mode emerges. 
As well-known, TA phonons possess a low frequency mode 
at the long wavelength, uj'^j^{q) cx which gives rise to 
the fluctuation of the FE polarization. Such polarization 
fluctuations are hybridized with the spin bosons and soften 
thus the transverse phonon behavior. 

Existing experimental data suggests that the exchange 
energy J falls off as a power law with the separation of 
the magnetic ions 

JiAn - rj) = JiMiRt + Ui) - {Rj + u,)|-^^^ (3) 

where 7 is in the range of 6 — 14 [16]. i?f is the bare 
value of the position of the atom at site i, and \R^ — Rj\ 
determines the lattice constant a (set here to 1). Uj is the 
displacements of site i. Generally, |ui| is small and does 
not destroy the lattice structure, (Juij/a ~ 10""^). We 
inspect the dominant term for J in the following two case. 

Longitudinal phonons. When the atoms are displaced 
along the chain one finds 

Ji.An - rj) = Ji,2[l - 71,2 • (u,^ - u|)] (4) 

with Cij being the unit vector connecting two sites i andj. 
A trillinear coupling between the phonon and spin is in- 
duced. One can easily check that (uf) = because of 
the local rotational symmetry of the spin-spin correlation 
(Si • Sj) in the helically ordered phase. Since uf is not 
involved in the DM interaction we do not consider it in 
the following discussion. 

Transverse phonons. For atomic displacements perpen- 
dicular to the chain, • 62 = we find for Ji,2(»'i — rj) 

JiMn - rj) « Ji,2[l - ^(u^ - uf)'] (5) 

which gives a TA phonon mode coupled to the spins with 

the bi-quadratic interaction —j1.2J1.2iu.j- — uj')^(Si • Sj). 
The effective spring constant for the TA phonon is kxA = 
jJ{Si ■ Sj). Furthermore, because of the negative Ji(Si • 
Sj) + J2{Si ■ S„j) the dynamical exchange-striction will 
harden the frequency of the transverse phonon mode as 

ijl=u,l{l + \F,\) (6) 
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where = k/M and Fs = [-71-^1(8^ • Sj) - 72J2(Sz • 
^rn)]/k- As the temperature decreases below the tran- 
sition temperature for magnetic order Tjv the spin-spin 
correlation function (S^ • S^) increases and so does the 
phonon frequency. Lowering further the temperature to 
the FE transition temperature TpE, an additional fre- 
quency hardening occurs due to the dynamical DM in- 
teraction [8,10]. The complete scenario is thus that the 
phonon frequency hardens at two onsets at Tjv and Tfe, 
a conclusion consistent with the experimental observation 
for Euo.75Yo.25Mn03 [12]. Assuming k ~ leV/A^ and 
JS'^ lOmeV [8] the frequency hardening can be esti- 
mated to be Slj/loo « 1%, which is in good agreement 
with the experimental data [12]. Phenomenologically, the 
exchange-striction suggests that in a Ginzburg-Landau 
(GL) theory for the coupling between the spin S and the 
transverse electric dipole P± terms of the form —aS'^Pj^ 
appear. As a consequence, P±_ and S condense at the 
same temperature due to the strong spin-lattice coupling 
a ^ J and the two transition temperatures merge, a con- 
clusion which is in line with the experimental observations 
in YMnOa [17]. There the electric dipole moment AP^, 
which is along the z direction obeys the same temperature 
dependence as the magnetic moment that is aligned in the 
ab{xy) plane with 120° structure below 80K. 



Electromagnon. — We split the atomic displace- 
ments into two parts: (i) the statical part Uj = {uq,0,0) 
driven by the DM interaction, and (ii) the dynamical part 
i5ui = {—Suf,duy,0) induced by the exchange-striction. 
As the softness of the system is due to the magnetic part 
we concentrate at first on the spin excitations. For the zx 
helical spins, it is convenient to rotate the spins locally (at 
each site) along its classical direction {Sf) 



Sf = Sf cos iQ + SfsiniQ 

sy = sf 

Sf = -Sf sin iQ + Sf cos iQ. 



(7) 
(8) 
(9) 



Disregarding the high-order terms of the interplay between 
the spins and the dynamical part of lattice displacements, 
i.e. using the standard linear-spin-wave approximation, 
we have 



A{q) = — Ji[cos(5 + -(1 + cosQ) cosg] 

- J2 [cos 2Q + + COS 2Q) cos 2q] 

A^S^sin^Q^^ 
+ —1 (2 — cosg) 



2k 



(11) 



B{q) 



^(1 — cosQ) cosg -h ^(1 — cos 2(5) cos2g' 



X^S^ sin " Q 
Ak 



■ cosq 



(12) 



H can be easily diagonalized by a Baguliubov transforma- 
tions. The energy dispersion of the spin-excitation reads 



(13) 



The effective spin anisotropy introduced by the spin- 
phonon(DM) interaction results in an energy gap of the 
spin-wave spectrum for non-coUincar spin ordering, i.e. 
^s{q = Q) ^ if Q ^ or TT. One can see in the fur- 
ther discussion that the spin fluctuation with the wave 
vector q K Q arc important in connection with the mag- 
netic softening of the transverse phonons. 

Now we turn our attention to the dynamical spin- 
phonon interaction. Retaining terms up to the second 
order in the quantum fluctuation, the spin-current model 
delivers the following coupling terms 

Hdm = -\S cos Q ^ (5<(5f+i - Sf) 

i 

-XS J2 ^u^Sf COS Qi+i - %i cos Qi) 

i 

= -XS cos Q^^6ugSg {cosq - 1) 

Q 

-XS^SuyS^^^Q{e^'Q -e'^'>^Q^)/2 (14) 

Sufj is hybridized with the spin at q±Q, but Su^ is coupled 
to S^ at q. As expected, SUg has the same long wavelength 
behavior as the magnons. No static displacement exists 
along the x direction, i.e. Su^ = 0. On the other hand, a 
uniform lattice deformation along the y direction {du^ ^ 
0) may occur due the hybridization between the electric 
polarization and the spin ordering [8] . After some algebra, 
we flnd for the polarization correlation functions 



H = Eo + J2^il)S^S+ + B{q){S-S^+S+Sf), 



M[a;4 - ^2(^2 + u;|) + a;2(a;2 - wf^)] 
1 

M[u,^-u,l + ^E,'=,±QGs{q')] 



where Eq = N[JiS'^ cos Q + J2S'^ cos 2Q 



(10) where LVp = sfw^'+U^^ is the frequency for the trans- 
verse phonon, J^^{q) = \l{A{q) - 2B{q)){}?S^ cos^ Q{\ - 

cosg))/fc']i/2^ and GM^Q) = {A{(i±Q)^2B{q±Q)){\- 
cos(g ± 2(5)/(a;2 — a;s(g ± Q)). At small wave vectors. 



, and 
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(/ ~ and lota{0) ~ 0, the TA phonon is decoupled from 
spins. Tlic antisymmetric DM interaction dominates over 
the spin-phonon coupling. Su^ is coupled via {Sq — Sq) 
to the rotation of the spin plane and the direction of the 
polarization along the chain. Additionally a uniform po- 
larization in the y direction is induced by the dynamical 
ME interaction and <C Su^^ISuq » possesses a low fre- 
quency behavior. The rotation mode around the z axis 
has Ll)'^ ~ V JSD if an easy-plane spin anisotropy D{Syy 
is introduced to the spin system. However, at a wave vec- 
tor close to the magnetic modulation vector, i.e. q ^ Q 
and ijjta{Q) 7^ both the symmetric and antisymmetric 
magnetoelectric interaction respond to the fluctuations of 
the polarization. Especially, in the direction parallel to 
the EE polarization P^, there is a low frequency range 
around = ^^(Q) where couples resonantly to light 
even \i D = 0. For finite D, assuming D ;» X'^S'^/k we 
observe nearly the same low-frequency behavior of the po- 
larization correlation functions u>^ « \/ JSD w . These 
conclusions are also qualitatively consistent with experi- 
ment observations (Fig. 8 in Ref. [13]). 

Spin- flip. — Recently, LiCu2 02 (5* = 1/2) has been 
found to be ferroelectric in the bc-spiral state at low tem- 
peratures [14]. In contrast to large spin multiferroics, in 
spin- 1/2 magnet the spin fluctuations may spontaneously 
reverses the local spin. According the spin-current model 
Eq.(2) the direction of the on site electric dipole mo- 
ment can also be completely reversed by the spin fluctua- 
tions. Large quantum fluctuations of the FE polarization 
^u-^ = — 2uq is induced by the hybridization between the 
phonon and the spin. 

Defining the vector of spin chirality as the average of the 
outer product of two adjacent spins Ci = (si x Si+i)/\si x 
in the RMnOa-typc multiferroics the direction of 
electric polarization is determined by the spin chirality. 
Reversing the direction of electric polarization does also 
reverse q. Clearly, the spin chirality has only two eigen- 
values, -|-1 and —1, and possesses long-Tange ferromagnetic 
order in the FE phase. Thus, Cj can be simply treated as 
the Pauli operator. According to the dynamical cxchange- 
striction Eq.(5) the interaction term involving the spin 
chirality has the structure —Jc{Q)ci ■ Ci^i- The x com- 
ponent of the spin-chirality operator acts as a direction 
reversal operator which can be traced back to the quantum 
fluctuation of the FE polarization. Considering the the 
dynamical DM interaction Eq.(14), the coupling term be- 
tween the spin and the spin-chirality in the spin-1/2 mul- 
tiferroics is given by J^i cf - sf) = J^i sf (cf-i - cf), 
which indicates that when the spin at site i is flipped, 
Si —^ —Si, the direction of spin-chirality Ci and Ci_i are 
also reversed, an observation consistent with the spin- 
current model and with the definition of the spin-chirality. 

For the one-dimensional spin-1/2 chain, the quantum 
model predicts a gapped spin-liquid state in the range 
of the frustration exchange parameters in LiCu202- The 
very existence of the magnetic helix state suggests that 



the quantum fluctuations is signiflcantly suppressed and 
the spins tend to recover a scmiclassical behavior. In the 
ground state of the spin system all spins point along their 
corresponding classical directions as in NaCu202, where a 
Ji — J2 spin model provides a good description of the helix 
state [15]. So the spin interaction can be simply given as 
— Js{Q)si ■ Sj where Q is taken as the pitch angle along 
the chain. 

Based on the above arguments an effective model that 
describes the interplay between the helical spin and spin- 
chirality has the form 

Hsc = - J^^-^'^i ■ + -^^^i ■ ^J') -^^^ii^i-i - Ci). (15) 

The Hilbert space can be considered as the tensor product 
space \i) Isf )s® jcf )c- The ground state of Hgc possesses 
the ferromagnetic order both for s and c, i.e. \g.s.) = 
\FM)s eg) \FM)c. Now let us consider the effect of the 
quantum fluctuation. If the spin at site i is flipped we 
have \si) = sflg.s.) = |sf)s \FM)c which is the ground 
state with the spin at site i being flipped. Noting that 
s-\g.s.) = |5f), ® \FM),, (sfflg.s.) = \g.s.), cf\g.s.) = 



\FM),(g>\cl) 
state ( hi 



{^i?\9.s.) 
we find 



Eo{q)+Es: 



\g.s.), if we apply Hgc to the 



cos^ 



sm.9 
- cos 6 



where the state \sq) {\cq)) is essentially a flipped 
spin (spin-chirality) delocalized across all the lattice. 
E^ = ^NJ,s^ - NJ.c^ + [JM + Jc(<z)]/2, J,{q) = 
2sJg{l — cosg), Jc{q) = 2cJc(l — cosg), and Egp = 
V[iJs{Q)-Mq))/2]^+l{p?, l{p) = 7(1 -cost?), cos 6 = 
{Jsip) - Jc{p))/'^Esp, and sin6' = ^{p)/Esp. Applying the 
|cg> \ ^ ( -sin 0/2 cos 0/2 \ / \sq) 
\sq) ) \ cos 61/2 sin 61/2 ) \ \cq) 
the Hamiltonian is brought in the diagonal form 



rotation 



Hsc = J2(^o{q) + Esp)\sq) + J2i^o{q) - E,p)\cg). 

q q 

Due to the spin-phonon coupling the spin and spin- 
chirality excitations are mixed. Two separated chan- 
nels are identified: the spin-channel \sq) and the phonon- 
channel jcg). In each channel, we have 



{s) + {£) = !. 



(16) 



Generally, the expected value of c is less then one due 
hybridization with spin excitations. A non-unitary c is the 
origin for a finite scattering intensity of polarized neutrons: 
For the cycloidal helimagnet, we have [18] 



(2) 



Ion - Ipff 
Ion + loff 



(17) 



where Ion{Iofj) is the reflection intensity of polarized neu- 
trons parallel(antiparallel) to the scattering vector. On 
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tlie basis of the experimental data for LiCu202 [14] wc 
infer (c) « 0.3. On the other hand, the magnitude of 
the ordered moment per magnetic copper site is 0.56yKB 
[15]. Together with the typical g-factor for Cu^+ in a 
square-planar geometry (g w 2) from Eq.(16) wc conclude 
(c) = 0.44 which is consistent with the previous estimated 
value. 

Summarizing, both the (symmetric) exchange-striction 
and (antisymmetric) DM interaction affect dynamically 
the magnetoelectric coupling in muitiferroics. At a small 
wave vector, the DM interaction determines the low- 
frequency behavior of the phonons. For a wave vector 
close to that of the magnetically modulated structure, the 
exchange striction induces fluctuations in the FE polar- 
ization, and additional low-lying mode parallel to the FE 
polarization emerges. For spin-1/2 muitiferroics, the effect 
of the quantum fluctuation is particularly large. The local 
polarization can be completely reversed by the spin fluctu- 
ation, and so does the direction of the on site spin-chirality. 
These findings are in line with experimental observations. 
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